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Abstract 

For n > d/2, the Sobolev (Bessel potential) space H n (R. d ,C) is known to 
be a Banach algebra with its standard norm || || n and the pointwise product; 
so, there is a best constant K n d such that ||/<7|| n ^ ^nd||/||n||<7||n f° r au 
/, g in this space. In this paper we derive upper and lower bounds for these 
constants, for any dimension d and any (possibly noninteger) n S (<i/2,+oo). 
Our analysis also includes the limit cases n — > (d/2) + and n — ► +oo, for 
which asymptotic formulas are presented. Both in these limit cases and for 
intermediate values of n, the lower bounds are fairly close to the upper bounds. 
Numerical tables are given for d = 1, 2, 3, 4, where the lower bounds are always 
between 75% and 88% of the upper bounds. 
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1 Introduction. 



The theory of Sobolev spaces contains a lot of inequalities which involve real con- 
stants; often, the classical arguments employed to prove these inequalities allow to 
infer the existence of such constants, but are unsuitable to evaluate them accurately. 
On the other hand, a precise knowledge of these constants is desirable for several 
reasons: apart from the intrinsic interest of the problem, there are many applications 
where a fully quantitative analysis relies on these numbers. 

The inequality analyzed in this paper refers to the pointwise multiplication in 
H n (R d , C) for any n > d/2. We are interested in the best constant K nd such 
that 

Wf9\\n^K nd \\f\\ n \\g\\ n 

for all f,g G H n (R, d , C), where || || n is the standard norm of this space (see Eq.s 
(1.2) (1.3) later on in this Introduction, and Eq. (2.1) in the next section). 
The constants K nd are relevant in relation to PDEs with polynomial nonlinearities, 
since they allow precise estimates on certain approximation methods and on blow up 
phenomena. To cite only one example, we refer to the semilinear heat equation in 
one space dimension discussed in [10]; here, an estimate on Ku has been employed 
to compute the error of the Galerkin approximate solutions, and the blow up times 
for certain initial data. 

Evaluating K nd for arbitrary n and d is a nontrivial task. For example, let the prob- 
lem be formulated in the variational language: maximize ||/<?|| n with the constraints 
1 1 /Ik = \\g\\n = 1; if n is integer one can write the corresponding Euler-Lagrange 
equations, but these form a cubic system of PDEs of order 2n for / and g. 
Due to the difficulty of the problem, one could be satisfied even if, in spite of the 
exact value of K nd , one had sufficiently close lower and upper bounds for it. Such 
bounds are proposed in this paper, for any integer d and (possibly noninteger) n e 
{d/2, +oo). Our upper bounds depend on an accurate use of the Fourier transform 
and of the convolution: the conclusion of this analysis is an inequality 

where K£ d is the sup on [0, +oo) of a function of hypergeometric type. This sup is 
easily evaluated, analytically in certain cases and numerically otherwise. 
The lower bounds we propose follow directly from the inequality that defines K nd , 
choosing for / and g appropriate trial functions: these often depend on one or two 
real parameters, so one gets the highest lower bound from the chosen functions max- 
imizing with respect to the parameters. In any case, this procedure gives inequalities 
of the form 

K nd ^ K nd , 

where K~ d depends on the trial functions: we will consider two specific choices, 
giving rise to what we call the "Bessel" or "Fourier" lower bounds. Both types 
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of bounds are expressible via special functions of hypergeometric type, or by one- 
dimensional integrals which are easily computed numerically. For given values of n 
and d, the best available estimate from below for K n d is obtained choosing for K~ d 
the highest between the Bessel and the Fourier bounds. For certain values of n only 
one kind of lower bound is easily computed, so one must be content with it. Our 
investigation also includes the limit cases n — > (d/2) + and n — > +00; the second 
limit requires the asymptotic analysis of certain integrals, which is performed via 
the Laplace method. To give an idea of our results, we anticipate some of them, 
i) For n -> (d/2) + , it is 



y/n - d/2 



where Md is an explicitly given constant (see the next section, Eq. (2.6)); on the 
other hand, denoting with K~ d a conveniently chosen Bessel lower bound, one finds 



K, 



nd 



3 - d/2 



~/ d s ~\ 
l + 0(n-- 



so, in this limit K+ d /K~ d -> ^2/1 > 0.816. 
ii) For n — > +00, it is 



Kd = T d 



(2/vsy 



n 



d/4 



n 



with Td another explicitly given constant (see Eq. (2.7)). On the other hand, 
denoting with K~ d an appropriate Fourier lower bound, one finds 



K, 



nd 



(5/3)^ (2/V 3)r 
-Id 



71/4 



71' 



d/4 



l + 0(-) 

n 



thus, K+ d /K- d - (5/3)^7-1/4 > 793 

iii) For d — 1,2,3,4 we have explored the whole interval n G (d/2, +00), choosing 
for each K~ d the most convenient Bessel or Fourier lower bound and comparing it 
with the upper bound K^ d ; for the sample values of n we have considered, K^ d /K~ d 
ranges between 0.750 and 0.880. A table of these upper and lower bounds is reported 
in the paper. 

iv) As previously said, K^ d is the sup of a hypergeometric-like function. Even though 
this is easily computed numerically, to avoid this burden one can use a majorant 
Kt/t ^ K\. We define K^f using only elementary functions of n; this bound 



reproduces correctly the asymptotic behavior of K^ d for n — > (d/2)' 
for 1 < d < 7 is very close to it on the whole range (d/2, +00). 



n 



+00, and 
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At the end of this Introduction we will give some details on the organization of 
the paper. Before speaking about this, we insert a few comments on some related 
literature. 

Connections with previous works. In our paper [8], we estimated the constants 
for more general inequalities related to multiplication in Sobolev spaces; in particu- 
lar, we discussed the constants K nad in the "tame" (or "Nash-Moser") inequality 

\\fg\\n < K nad m&x(\\f\\ 

for d/2 < a < n and /, g e H n (R d , C); here || || a is the norm of H a (R d , C). (The 
cited work is partly related to the previous one [7], and to the subsequent one [9] 
on the tame functional calculus in Sobolev spaces). In the special case n = a, the 
inequality written above coincides with the inequality of the present paper. 
For arbitrary d, a, n, in [8] we derived upper and lower bounds for K nad . The lower 
bounds were of the Bessel and Fourier types also considered here (with no analysis 
of the limit n — > (d/2) + , and a discussion of the limit a fixed, n — > +oo, of course 
different from the present limit n — > +oo; some explicit formulas of [8] for these 
lower bounds are replaced here with equivalent, but simpler versions, and we also 
give some new formula). 

The upper bounds for K nad were obtained by a different method than the present 
one for K nd \ furthermore, if the upper estimates of [8] are applied with n = a they 
are found to be rougher than the present ones on K nd . 

The method we use here to get the upper bounds refines an idea which appeared in 
[13] in relation to the multiplication in the space H n (T, C), where T := R/ (2ttZ) is 
the one-dimensional torus. The author of [13] was not interested in a precise estimate 
of the constant for multiplication, so he inserted in his argument some majorization 
which, although unnecessary, simplified the proof of the convergence of a series; the 
upper bound on the constant for the multiplication in H n (T, C) arising from this 
simplification behaves like const. x2 n for large n (see page 294 of the cited paper). 
Here we replace the one- dimensional torus with R d , and the Fourier series with the 
ci-dimensional Fourier transform. The literal translation of the technique of [13] 
in our framework would give again an upper bound for K nd behaving like 2 n for 
n — > +oo; on the contrary, here we use only the strictly necessary majorizations 
and finally obtain the bound K^ d involving a hypergeometric function, which as 
explained behaves like (2/v^3) n n~ d//4 for n — > +oo and is accurate for small n as 
well. 

Organization of the paper. In Section 2 we state precisely all the results about 
the previously mentioned upper and lower bounds for K nd . As a preparation for 
the proofs, in Section 3 we write a list of known identities frequently cited in the 
sequel, on the following subjects: radial integrals, radial Fourier transforms, hyper- 
geometric functions, integrals with three Bessel functions and the asymptotics of 
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Laplace integrals (the last two topics are also treated in the Appendices A and B). 
In Section 4 we prove all statements about the upper bounds K+ d . In Sections 5 and 
6 we prove all the results about the Bessel and Fourier lower bounds, respectively. 
In the remaining part of this Introduction, we fix some notations and definitions 
employed as standards throughout the paper. 

Basic notations on R d and Fourier transforms. We consider an arbitrary 
space dimension d; the running variable in H d is x = (x±, ...,#<*), and k = (ki, k d ) 
when R d is interpreted as the "wave vector" space of the Fourier transform. We 
write • and | | for the inner product and the Euclidean norm of R d (so that | x\ = 
\/xi 2 + ... + x d 2 , \k\ = \/k\ + ... + k d 2 , k • x — kiXi + ... + k d x d ). 
We denote with J 7 , T~ x : S'(R d , C) — > S'(R d , C) the Fourier transform of tempered 
distributions and its inverse, choosing normalizations so that (for / in L 1 (R d , C) ) it 
is Ff(h) = (27r)- rf / 2 J Rd dx e~ ik * x f(x). The restriction of T to L 2 (R d , C), with the 
standard inner product and the associated norm || \\ L 2, is a Hilbertian isomorphism. 

Sobolev spaces. For real n ^ 0, let us introduce the operators 

S'(R d , C) - S'(R d , C) , g » VT^A n g := T' 1 (a/TT^ (1.1) 

where sj\ + |k| 2 ™ means the function k E H d h- > \J\ + \k\ 2 ™. The n-th order 
Sobolev (or Bessel potential [3]) space of L 2 type and its norm are 

H n (R d , C) := {/ G S'(R d , C) VT^A n f e L 2 (R d ,C)} = (1.2) 

= {feS\R d ,C) | v^Tpf V/GL 2 (R d ,C)} , 

\\f\\ n := ||v / T^A " /|| i2 = || " ^/ • (1.3) 

For n integer, these definitions imply 

H n (R d , C) = {/ e S'(R d , C) | V m / G L 2 (R d , ® m C d ) Vm G {0, n} } (1.4) 

where 

V m / := (9 Al ,..., Am /) (Al ,..., Am)e{1 ,..., 4m (1.5) 

and (9 Ai is the distributional derivative with respect to the coordinate x\ t . The 
statement V rn f G L 2 (R d , ® m C d ) means that 

> > ; / dx\(d Xl ,...,u)(x)\ 2 --=\\v m f\\ 2 L2 , (i.6) 



Ai ,...,A m — 



and the norm (1.3) can be written as 



m=0 



E ( I ) W^fWh ■ (1-7) 



Other notations. Some useful functions. The Pochhammer symbol of a E R, 
£ G N, is 

(a) e :=a(a + i)...(a + e-i) . (1.8) 
The semifactorial of an odd mGNis 

m\\ := 1.3....(m-2)m , (1.9) 

and we also intend (—1)!! := 1. We refer to [1] [5] [14] as our standards for special 
functions. In this paper, we frequently use the Gamma function and its logarithmic 
derivative ip(w) := T'(w)/T(w); for future reference, we write here their properties 
more frequently employed in the sequel. These are: the shift formulas 

F(w + 1) = wT(w) , (1.10) 

il>(w + 1) = + - ; (i-ii) 

w 

the special values 

r(l/2) = v^F , r(l) = 1 , V(V2) = -7. - 2 log 2, ^(1) = - 1e (1.12) 
(with 7 B the Euler-Mascheroni constant); the duplication formula 

o2wj-1 

r(2w) = — =- r(«; + i/2)r» ; (1.13) 

V 71 " 



the identity 



/ 

Jo 



= f„ r7 >,>0. (1.14) 

(i + up r( 7 ) 



Another function of which we make wide use is the Gaussian hypergeometric function 
2F±(a, (3, 7; w) = F(a, (3,~f;w). We are especially interested in the function 

F n d '■ [0, +00) — >■ (0, +00) , u F„ d (w) := F ( 2n - ^,n,n + ^; -u J , (1.15) 
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rf€N\{0}, nG(rf/2,+oo). 
This function has the equivalent representation 

F nd (u) = - — l —— F ( n, - + - - n,n + -; — ^— ) , (1.16) 

following from a familiar Kummer transformation (see Sect. 3, where we return to 
some statements appearing here); we also mention the special case 

nd[) (n + l/2) e £\ (1 + «)»+< 2 2G 1 j 
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As it often happens dealing with Sobolev spaces, a central role in our considerations 
is played by the functions 

G nd :R d ^C, k^G nd (k) := (1+ | fc | 2)n i (1-18) 

Qnd '■ R d — ► C , g nd := T 1 G nd ■ (1-19) 
It is clear that g nd G H n (R d , C) if n > d/2; explicitly, one has [3] [6] 

9n d (x) = 2 '„' lr(n) K n _ d/2 (\x\) (1.20) 

for x £ R d ; here K u are the modified Bessel functions of the second kind, or Mac- 
donald functions. 



2 Description of the main results. 

Let d e N\ {0}. For (integer or noninteger) n > d/2, the space H n (R d , C) is known 
to be a Banach algebra under the pointwise multiplication: see, e.g., [2]. 

2.1 Definition. For n > d/2, we put 

K nd : = min { K > | ||^|| n < K\\f\\ n \\g\\ n for all f,ge H n (R d , C) } (2.1) 

and refer to this as the best (or sharp) constant for the multiplication in H n (R d , O). 
In the sequel we present our upper and lower bounds on K nd . 

Upper bounds on K nd . These are given by the following proposition, to be proved 
in Sect. 4. 

2.2 Proposition, i) For all n > d/2, it is 



K nd ^K+ d := sup y nd {u) . (2.2) 
Y «e[o,+oo) 

y nd : [0, +oo) -> (0, +oo) , y nd {u) ■= ^/TnanV 1 + ^ Fnd{u) ' (2 ' 3) 

with F nd as in Eq. (1.15) or (1.16). y nd is bounded, and its boundary values for 
u = , u — > +oo are 

r(2n-d/2) r(n + l-d/2) 

" (47r)<V 2 r(2n) ' ^ +00 J " 2<*-%<*/ 2 (n - d/2)r(n) ' 1 J 
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ii) For d/2 < n ^ d/2 + 1/2 the function y n<d is increasing, so that 

1 



K 



nd 



OO 



T(n + 1 -d/2) 



2 d/2-i/2 n d/4 y ( n _d/2)r(n) ' 

For fixed d and n — > (d/2) + ; this implies 



(2.5) 



1 + 0(n 



^/ra - d/2 
ra) For /raed d and n — > +oo ; it is 



d 



2<*/2-i/2 7r d/4 v /r(d/2) 



(2.6) 



l + OR 



n 



(2/v/3) r 



d/4 



l + OR 



n 



gd/4+1/4 

2 d ir d / 4 ' 



(2.7) 



Of course, in Eq.s (2.4) and (2.5) we could write F(n + 1 - d/2)/(n - d/2) = 
T(n — d/2); the expression in the left hand side has been preferred to handle the 
limit n — > (d/2) + . Similar choices have been made for other formulas in the sequel. 

"Bessel" lower bounds on K n d- The general method to obtain lower bounds on 
this constant is based on the obvious inequality 



K nd ^ 



Wfgl 



\9\ 



(21 



for all nonzero f,g e H n (R d , C); this gives a lower bound for any pair of "trial 
functions" f,g. Inspired by [8], we choose for / and g the function 



g\nd{x) ■■= gnd(^x) 



(2.9) 



where A G (0, +oo) is a parameter and g nd is defined by Eq. (1.19). By comparison 
with that equation, we find 



g\nd — F G\nd , G\ n d(k) 



1 



A d (l + \k\ 2 /\ 2 ) n ' 



(2.10) 



To give a lower bound for K nd in terms of these functions simply amounts to compute 
llfl'ArKflU) llfi'AndlU- These norms were already calculated in [8]; here we give them in 
a more simple and complete form, and add an analysis of the limit case when n is 
close to d/2. All these facts are described by the forthcoming proposition, to be 
proved in Sect. 5. 



hi 



nd\\n 



2.3 Proposition, i) For all n > d/2 and A > 0, it is 

7 



\9\nd\ 



(2.11) 



whence 

K nd > K* d := sup JCW • (2-12) 

A>0 

The norms in Eq. (2.11) are given by 

ll2 n d / 2 T{n + 1 - d/2) d , x2 , , rtin , 



7T 



d/2 



" r(d/2)r(2n)A" x (2 - 14) 



wii/i F n d as in Eq.s (1.15-1.16); 



2 2 = 7r rf / 2 r 2 (2n-d/2) " "^ 1/2 (n), (d/2 + 1/2 - n), (n) ro (d/2 + 1/2 - n) 



m x 



r(d/2)P(2n)A d (n + l/2)e£\ (ra + l/2) m m! 

r(d/2 + £ + m)r(n-d/2) d . . , A ^ 

x w i o\ \ F -n,- + £ + m,n + £ + m;l-4A 2 

/or n — | — | integer . (2-16) 
Lei d/2 < n < d/2 + 1/2. JTiera, /or all \ > it is 

so i/ia£ 

^(A) > ^ B (A) := , (2.18) 

tf£ > K% := sup J#j fl (A) . (2.19) 

A>0 



Here: 



f (X)- ^ /2r2 ( 2 ^/ 2 ) \ P 2 r(n+l-d/2) d 

/UA) - L nd — fw — ( ' 2' n ' 1_4A )+ (2 - 20) 

„ ^ T(2n+1 — d) d d ,, 9 . 

" P -«- r(2n-W F{ ~"' 2 ' 2 " " 2 ; 1 " 4A )+ 
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I> + l/2)r(n + l-d/2) 
Fnd - ^Y{2n-d/2) ' [Z - ZL) 

r(n + i/2)r(d/2 + i- 7i) r g nd */ p nd > g nd; 

^ nd ■ r(n)r(l/2 + d/2 - n) ' 9nd ' \ - (n - d/2) i/ P nd < Q nd . 

(In the above definition of Q n d one should intend T(0) := oo ; so that Q nd = for 
n = d/2 + 1/2). For any fixed d, A and for n — > (d/2) + , it is 



^( A ) = \/^^=^ l + , (2-22) 



2 Afw 



3 - d/2 



2' 



with M d as in the asymptotic expression (2.6) for the upper bound K^ d (note that 
y^2/3> 0.816/. 

As clarified in the sequel, the Bessel lower bounds are less interesting for large n; 
therefore, it is not worth to determine their asymptotics for n — > +oo. 

"Fourier" lower bounds on K n( i- Another choice for the trial functions amounts 
to choose for / and g the function 

fp*d(x) := e** 1 e -W 2 )W a (2.23) 

where the "Fourier character" x — > e' ipxi is regularized at infinity by a Gaussian 
factor (we take this hint from [8], but we develop it in a different way). 
As we will see, this choice is especially interesting for large n. The Sobolev norm 
of any order n of this function can be expressed using the modified Bessel function 
of the first kind I u , the Pochhammer symbol (1.8) and the semifactorial (1.9). Our 
results on the Fourier lower bounds are contained in the forthcoming proposition, 
to be proved in Sect. 6. 

2.4 Proposition, i) Let n > d/2. For allp,a > 0, it is 
hence 

K nd > K F nd := sup JC(P, °) ■ (2-25) 

p,cr>0 

For all p,a > 0, it is 

WM\l = J^-, J^dpp^il+pre-^h^p) ; (2.26) 
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in particular, for n integer it is 



£ 

n 



ii/^-^eee(';j(;j(1)^x 

1=0 j=0 g=0 V / V J / V » / 



x (d/2 - l/2)^.p 2j - 2 V^- J - d / 2 . (2.27) 
ii) Fix the attention on the "special" lower bound 

KU~XZ{p=^, =$; (2.28) 

then 

'FF (5/3) 1 / 2 (2/VS) B L 



7^2<_f V-/ "/ rp 

-^nd - -71/4 J rf 



7 l/4 ^ n d/4 



n J 



l + 0(-) /orn^+oo, (2.29) 



with T d as in the asympotic formula (2.7) for the upper bound (note that (5/3) ^ /7 1 / 4 
> 0.793;. 

Remark. The result (2.29) depends on the asymptotic analysis of a Laplace inte- 
gral. The values for (p, a) in Eq. (2.28) have been chosen because they simplify this 
analysis, and give rise to the term (2/\/3) n n~ d ^ 4: also appearing in the asymptotics 
(2.7) for the upper bound. One could discuss the asymptotics of J^(p,a = c/n) 
for arbitrary choices of p and c in (0, +00); however, this generalization complicates 
the implementation of the Laplace method and, in comparison with (2.29), yields 
no sensible increase of the dominant term. 

Table of the upper and lower bounds on K nd for d = 1, 2,3, 4 and some 
test values of n. This is Table 1, which has been constructed using the upper 
bounds K^ d given by Prop. 2.2, and choosing conveniently one of the lower bounds 
K^ d) K^f , K^ d) K^f in Propositions 2.3, 2.4; the chosen lower bound is generally 
indicated with K~ d , and its type is specified within the table. We have chosen the 
values of n within a very wide range, from d/2 + 10~ 4 to d/2 + 120; for a better 
appreciation of the discrepancy between the upper and lower bounds, instead of K~ d 
we have reported the ratio K~J K^ d . 

To compute K+ d , we must find the sup of the function S^ nd in Prop. 2.2, which is 
given explicitly by item ii) of the same proposition for d/2 < n ^ d/2 + 1/2, and 
must be computed directly from the function S^ nd in the other cases; we have done 
this numerically in most cases, and sometimes analytically: some examples are given 
in Sect. 4. For large n, the numerical search for the maximum of 5^ nd has been done 
starting from u — 1/2, as suggested by item iii) of Prop 2.2. 

Concerning K~ d , we have always chosen for it the most convenient between the lower 
bounds in Propositions 2.3 and 2.4 (i.e., the highest one or, in some limit cases, the 
most easily computable). 
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As for the Bessel lower bounds, for n sufficiently distant from d/2 we have computed 
numerically the function A — > J£^(A) anc ^ ^ s maximum K^ d . For n very close to 
d/2, this computation is very difficult because the integrals in converge too 
slowly; in this case, we have turned the attention to the function A — > ^J B (A) and 
estimated numerically its maximum K®f . 

Concerning the Fourier lower bounds, for n not very large we have determined K^ d 
maximizing numerically the function (p,cr) — > J^(p, a); for very large n, we have 
turned the attention to the bound K^f which is easily computed numerically. 
The Bessel lower bounds are generally higher than the Fourier ones for small n; the 
contrary happens for large n. 

A more accurate n — > (d/2) + asymptotics for K^ d . This is introduced for 
the reasons explained in the next paragraph. For the sake of brevity, let us put 

n d :=n-^; (2.30) 

in place of Eq. (2.6), we propose a higher order expansion 

^l-JWOCnJ)!, N d : = tm±J± . (2 .31) 



This is derived from the explicit expression (2.5) of K* d , inserting therein the ex- 
pansions 

T(l + n d ) = T(l) + r'(l) n d + 0(n 2 d ) = l- lE n d + 0(n 2 d ) , (2.32) 



r(n) = T(^) + r'(^) n d + 0(n d ) = T(^) [l + ^) n d + 0{nl)] , 
(recall that Y'(w) = F(w)if>(w), and use Eq. (1.12)). 

"Elementary" upper bounds K^J~. The results (2.31) (2.7) on the asymptotics 
of K+ d in the limits n — > (d/2) + , n — > +oo suggest a way to build new majorants 

Kd > Kd >K nd , ne (d/2, +oo) , (2.33) 

that are presented hereafter. Even though less precise than the + upper bounds, 
the ++ bounds have the advantage of being elementary functions of n; we will show 
that they are very close to the + bounds on the whole interval (d/2,+oo) up to 
d = 7, and fairly close to them up to d = 10. For any d, the elementary ++ bounds 
reproduce the asymptotics (2.31) (2.7) of the + bounds at the leading order. 
In order to construct K^ d , we first define a function n G (d/2, +oo) h- > z nd through 
the equation 



,d. 



, M± u _ v±f (1 + Vdnd) + Td nnf + Znd - 

8 Jn d \ n / \ n / n J 



'nd n rf/4 
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V d := log(^) + \ + l~ N d, n d as in (2.30) . (2.34) 

This equation is easily solved for z nd . From the explicit expression for z nd and from 
the asymptotics (2.31) (2.7), one gets 

Znd = 0{^/n~ d ) for n -> (d/2) + , z nd = 0(1) for n -> +oo ; (2.35) 

the coefficient V d is defined as above just in order to give the first one of these 
relations. 

On account of Eq.s (2.35), for fixed d the function n — > z nd is bounded on the 
interval (d/2, +oo); this ensures the finiteness of 

Z d := sup z nd . (2.36) 

ne(d/2,+oo) 

Now, putting 



■ 3/2 . ~ ^~ 

77, J 

we see from (2.34) that K+ d ^ K+ d + . From Eq.s (2.37) and (2.31) (2.7), we also infer 



K ++ K ++ 1 

— ^- = l + 0(n d ) forn^(rf/2)+, _»sL = i + 0(_) for n -> +00 . (2.38) 

The forthcoming Table 2 reports, for 1 ^ d ^ 10, the numerical values of the 
constants Z d in Eq. (2.36) and of the quantities 

& d := sup (2.39) 

ne(d/2,+oo) 

The table has been constructed in this way. First of all, for each d in the above 
range the function n 1— > £ nd defined by (2.34) has been plotted (expressing z nd in 
terms of K+ d and evaluating the latter numerically); from the graph of n 1— > z nd , the 
sup Z d has been evaluated. Secondly, for the same values of d the ratio K^ d / ' K* d 
has been plotted as a function of n, and its sup d has been evaluated from the 
graph. 
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Table 1. Bounds K~ d < K nd < K+ d for d = 1, 2, 3, 4 and n - d/2 = 10~ 4 , 
lO -2 , lO -1 , 1/4, 1/2,1,3/2,3,6, 15,30, 60, 120. (The symbol stands for 
one of the types bb, b, f, ff, indicated below.) 



d = 1 



n 




i+io- 2 


i+io- 


3/4 


1 


3/2 


2 


7/2 


13/2 


31/2 


61/2 


121/2 


241/2 


K d 


56.5 


5.69 


1.90 


1.30 


1.00 


0.852 


0.814 


0.834 


1.07 


3.09 


22.4 


1410 


6.63 xlO 6 


< d 


0.816 

(BB) 


0.818 

(BB) 


0.824 

(BB) 


0.834 

(B) 


0.842 

(B) 


0.810 

(B) 


0.777 

(B) 


0.766 

(F) 


0.787 

(F) 


0.794 

(F) 


0.794 

(F) 


0.789 

(FF) 


0.791 

(FF) 



d = 2 



n 


l+io- 4 


l+io~ 2 


l+io- 1 


5/4 


3/2 


2 


5/2 


4 


7 


16 


31 


61 


121 


K d 


39.9 


3.99 


1.27 


0.798 


0.565 


0.428 


0.378 


0.332 


0.361 


0.831 


5.08 


269 


1.07xl0 6 


K~ d 

K d 


0.816 

(BB) 


0.817 

(BB) 


0.826 
(BB) 


0.844 

(B) 


0.865 

(B) 


0.842 

(B) 


0.811 

(B) 


0.752 

(F) 


0.772 

(F) 


0.788 

(F) 


0.794 

(F) 


0.786 

(FF) 


0.789 

(FF) 



d = 3 



n 


f+io- 4 


f+io- 


f+io- 1 


7/4 


2 


5/2 


3 


9/2 


15/2 


33/2 


63/2 


123/2 


243/2 


K d 


22.6 


2.25 


0.692 


0.421 


0.283 


0.198 


0.164 


0.128 


0.120 


0.223 


1.15 


51.2 


1.71 xlO 5 


K~ d 

K d 


0.816 

(BB) 


0.817 

(BB) 


0.826 

(BB) 


0.847 

(B) 


0.875 

(B) 


0.858 
(B) 


0.830 

(B) 


0.763 

(B) 


0.759 

(F) 


0.781 

(F) 


0.788 

(F) 


0.782 

(FF) 


0.787 

(FF) 



d = 4 



n 


2+io- 4 


2+io- 2 


2+10- 1 


9/4 


5/2 


3 


7/2 


5 


8 


17 


32 


62 


122 


K d 


11.3 


1.12 


0.340 


0.202 


0.130 


0.0857 


0.0678 


0.0473 


0.0389 


0.0590 


0.259 


9.72 


2.73 xlO 4 


K~ d 

K d 


0.816 

(BB) 


0.817 

(BB) 


0.826 
(BB) 


0.849 

(B) 


0.880 

(B) 


0.867 

(B) 


0.842 

(B) 


0.779 

(B) 


0.750 

(F) 


0.775 

(F) 


0.785 

(F) 


0.778 

(FF) 


0.785 
(FF) 



Table 2. Constants and (for the elementary upper bounds K^ d ). 



d 


i 


2 


3 


4 


5 


6 


7 


8 


9 


10 


z d 





0.00925 


0.0458 


0.0782 


0.105 


0.122 


0.128 


0.125 


0.115 


0.102 


e d 


1.041 


1.039 


1.044 


1.044 


1.044 


1.044 


1.049 


1.105 


1.197 


1.363 
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3 Some background. 



In this section we review some known facts, frequently cited in the rest of the paper 
to prove the statements of Sect. 2. 

Some d-dimensional integrals. We frequently need to compute integrals of func- 
tions on H d which depend only on the radius | | (radially symmetric functions), or 
on the radius and one angle. In this case, we use the formulas 



L 



r 2 7r ' 

^ d dx M) = r#j I drrd "V(r); (3.1) 

27r d/2-l/2 

*°x(\x\,Vx) = m2 _ 1/2) * (3-2) 



x 



/» + 00 /*7T 

/ drr d ~ l d9 sm6 d - 2 X (r,rcos9) (d > 2; rj e R d , \rj\ = 1) , 
Jo Jo 



holding for all (sufficiently regular) complex valued functions ip on (0, +oo) and \ on 
(0, +oo) x (0,7r). (When writing the analogous formulas for integrals on the "wave 
vector" space (R. d ,dk), the radius r will be renamed p). 

Radial Fourier transforms. Consider two (sufficiently regular) radially symmet- 
ric functions 

/:R*-C, x^f(x)=<p(\x\), F:R d ^C, k - F(k) = $(|fc|) ; (3.3) 

the Fourier and inverse Fourier transforms J-f, T~ l F are also radially symmetric, 
and given by [4] 

i /-+oo 

W)(*) = ]fcP^=T I drr d / 2 J d/2 ^(\k\rMr) , (3.4) 

^- l F){x) = - WI - 1 / rfpp rf / 2 J d/2 _ 1 (|x|p)<l.(p), (3.5) 
Fl -'o 

where are the Bessel functions of the first kind. As anticipated, the latter formula 
allows to infer Eq. (1.20) of the Introduction; in this case, Eq. (3.5) is applied with 
$(p) = 1/(1 + p 2 ) n and the corresponding integral over p is given in [14], page 434. 

Hypergeometric function. As anticipated, in this paper we use extensively the 
function F(a, /3, 7; w); we always interested in real values of the parameters a, (3, 7 
and of the argument w. For future citation, we report here some properties of F. 
First of all, we cite: the symmetry property 

F(a, (3, 7; w) = F((3, a, 7; w) ; (3.6) 
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the special values 

F(a,/3, 7 ;0) = 1 , (3.7) 
F(a, (3, r , 1) = ^ 7)r(7 ~ a ~^ for 7 > a + (3, 7 ^ 0, -1, -2, ... ; 

r(7 - a)r(7 - p) 

the particular cases 

F(a,[3,f3;w) = (l-w)- a , (3.8) 



(^M - m )£ ^ 

F(a, — m, 7; u») = > — lor m e N 

^— ' (7j£ £! 



^=0 



(3.9) 



Secondly, we recall that 

F(a, (3, 7; w) = f7^Y^~Tgj ds^-\l- s)^-\l-ws)- a > (3.10) 

for 7 > /3 > 0, w < 1 , 

F(a,/3, 7 ;l-w)= ^_ ^ jf duu^l + ^(l + «;«)""> (3.11) 

for 7 > /3 > 0, w > 

((3.11) follows from (3.10) with a change of variable s — u/(l + u)). 
Thirdly, we mention the differentiation formula 

■^-F(a, (3, 7; w) — — F(a + 1, /3 + 1, 7 + 1; «;) ; (3.12) 
aw 7 

this formula, combined with the positivity statement in (3.10) implies 

4-H<x, P, T, w) > for a > 0, 7 > (3 > 0, iu < 1. (3.13) 
aw 

Finally, we report the Kummer transformations 

1 w 

F(a,(3,r,w) = - _F(/3, 7 -a, 7 ; -) , (3.14) 

(1 — w)p w — 1 

F(a, (3, 7; w) = (l- wy-^F^ - a, 7 - /3, 7; w) ; (3.15) 

the first one allows to pass from the form (1.15) to the form (1.16) for F nd . The 
positivity of F nd is granted by (3.10). The expression (1.17) of F nd for n — d/2 — 1/2 
integer follows from (1.16) and (3.9). 

An integral involving Bessel functions. In Sect. 4 we will use the integral 

r+00 

I, v {h):= / drr^ +1 J^hr)Kl /2 {r) (// > -1, v > 0, h > 0) , (3.16) 
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involving a Bessel function of the first kind J M and the square of a Macdonald 
function K u / 2 . It is 

2^(^ + ^/2 + 3/2) h X ^- i7j 

xF(/i + 1/ + 1, n + z//2 + 1, n + z//2 + 3/2; -h 2 /A) . 

This result is probably known, but it is not easy to trace it in the most common 
tables on integrals of Bessel functions; for this reason, the proof of (3.17) is given in 
the Appendix A. 

Laplace integrals. The classical theory of these integrals is widely employed in 
this paper, to discuss the n — > +oo asymptotics of our bounds on K n d- 
By a standard Laplace integral, we mean an integral depending on a parameter n, 
of the form 

,6 

L(n) : = / dt tf(t) e~ n ^ , (3.18) 
Jo 

under the following assumptions: 

< n < n < +oo , < b < +oo ; (3.19) 
(p £ C 1 ((0, 6), R) , <p'(t)>0 Vte(0,6), hm + v?(t) = 0, 

r b 

d e C((0, 6), R) , / dt |tf(t)| e- nv(<) < +oo for all n as above . 

Jo 

Here and in the sequel, ' is the derivative; we shall also put 

£:=4eC7((0,&),R) . (3.20) 

The Laplace method gives the n — > +oo asymptotics of L(n), using the idea that 
the major contributions to this integral should come from the regions close to the 
minimum point of ip, i.e., to t = 0. (In certain cases, this asymptotics gives a fairly 
good approximation of L{n) also for non large values of n). The asymptotic behavior 
of L(n) is described by the following proposition (see, e.g., [11]; for uniformity of 
language, the proof is reviewed in the Appendix B). 

3.1 Proposition. Suppose that conditions (3.19) hold, and that 

£(t) = PMt)^- 1 + O^tr- 1 ) for t - + , (3.21) 
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where £ G {1, 2, ...} ; P ± , P t _ x G R, < a ± < a 2 < ... < a e . Then 

L W = E P ^ + °(i) fom^+oc. (3.22) 



=o 



More on Laplace integrals. By a general Laplace integral, we mean an integral 
depending on a parameter n of the form 

A(n) := f ds 0(s)e- n * (s) , (3.23) 

J a 

where 

< n < n < +oo , -oo < a < c < +oo , $ G C 1 ((a, c), R) , (3.24) 
9 G C((a,c),R), / rfs |0(s)|e- n$(s) < +oo for all n as above . 

J a 

Under suitable conditions on $, A(n) can be expressed in terms of one or more 
standard Laplace integrals. As a first example, suppose 

a > -oo , $'(s) > for all s G (a, c) , $(a) := lim $(s) > -oo (3.25) 

(the limit certainly exists by the monotonicity of $, but it could be — oo); then 

A(n) = e~ n<I>(a) L(n) , (3.26) 

L(n) as in (3.18) with 6 := c - a, p(t) := $(a + t) - $(a), tf(i) := 0(a + i) . 
Similarly, if 

c<+oo, $'(s) < for all s G (a, c) , $(c) := lim $(s) > -oo , (3.27) 
we can write 

A(n) = e- n * (c) L(n) , (3.28) 
L(n) as in (3.18) with b:=c-a, tp(t) := $(c - i) - $(c), tf(i) := 0(c - t) . 
As a final example, suppose 

| forsg/i (/iG(a,c)); (3.29) 

then we can write 

A(n) = e-" $w [L-(n) + L + (n)] , (3.30) 

L T (n):= dt^(t)e- n<p *V , 6" := /i - a, b + :=c-h, 
Jo 
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(p^(t) :=$(h^t) , ^(i) := 6(hTt) for i e (0, &*) , 

and L T (n) are standard Laplace integrals. 

In all the previous examples, after reexpressing A(n) in terms of standard Laplace 
integrals one should expand in powers of tp or (p T the functions C, '■= 6/<p' or 
£ T := 9 T /ip T '. Assuming sufficient smoothness for and $, the coefficients of 
these expansions can be expressed directly in terms of the derivatives of and $ 
at s = a, c or h, respectively [11]. In the third example, where $ has its minimum 
at an inner point h of (a,c), there is typically an alternation of equal and opposite 
coefficients in the expansions of £~ and £ + ; this yields some cancellation effects in 
the expansion of L~(n) + L + {n). 

4 Proofs for the upper bounds on K n d- 

Let us write F * G for the convolution of two (sufficiently regular) complex functions 
F, G on R d , given by 

(F*G)(k):= f dh F(k - h)G(h) . (4.1) 

We have 

Hf9) = j^ ir2 Ff*F9 (4.2) 

for all sufficiently regular functions / and g on H d (and in particular, for /, g as in 
the forthcoming Lemma). 

4.1 Lemma. For all n > d/2, it is 



K nd < /sup S nd (k) , (4.3) 

Y k£R d 

where 

S nd (k) := (1 (G nd * G nd ) (k) (4.4) 

and G nd (k) := 1/(1 + lA;] 2 )" for all k E R d , as m Eq. (1.18). 
Proof. Consider any two functions /, g e H n (R d , C). Then 

n/c= ^/Mi+^i^)wi 2 = ^/ Rd ^(i+^ 2 riw*^)wi 2 - (4.5) 

On the other hand, by making explicit the convolution we find 

(Ff*Fg)(k) = [ dhTf{k - h)Tg{h) = (4.6) 



18 



= / dh J" „ (Vl + I* - ^l 2 " J7(* - h) Vl + W Tg{h)) 

Jn d y/l + \k- h\ 2 y/1 + \h\ 2 v / 

Now, Holder's inequality \ J dh U(h)V(h)\ 2 ^ dh\U(h)\ 2 ^(j dh \V{h)\ 2 ^j gives 

\(Ff*Fg)(k)\ 2 ^C nd (k)P(k) , (4.7) 
f dh 

Cnd{k) '= L {i + \k-h\ 2 Y{i + \h\ 2 Y = {Gnd * Gnd) {k) ' 

P(k) := f dh(l + \k- h\ 2 ) n \Tf{k - h)\\l + \h\ 2 Y\Fg{h)\ 2 . 
Jn d 

Inserting (4.7) into Eq. (4.5) we get 

11/9111 ^wrL dk{l + \ k ^ n °^ p ^ ^ ( 4 - 8 ) 

<- ( sup ^V k l )n C nd (k)) [ dkP(k) = ( sup S nd (k)) [ dkP(k) . 

But 



/ dkP(k) = (4.9) 

= dk(i + \k\ 2 r\Ff(k)\ 2 ) ( dh(i + |/f n^(/i)| 2 ) = ii/n 2 ||^| 2 t , 



so we are led to the thesis. o 
4.2 Lemma. For n > d/2 and k G R d , it is 

s nd {k) = yj^-) , (4.io) 



4 

where 5^ nd is the function in Eq. (2.3) of Prop. 2.2. 

Proof. Let us recall that G nd is the Fourier transform of the function g nd , already 
considered in Eq.s (1.19) (1.20). We have 

S nd {k) = 1 J 2 ly {Tg nd * Tg nd ) (k) = 1 ^ )d \J (fg 2 nd ) (k) . (4.11) 

But g 2 d is a radially symmetric function, whose explicit expression in terms of the 
Macdonald function is given by (1.20). We insert this expression in the formula 
(3.4) for the radially symmetric Fourier transform and obtain 

i r+oo 

(*) = 2 2 ^ 2 T 2 (n)\k\^ J, dr r2n ~ d,2j ^\ k \ r ) K n-^) ; (4.12) 
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the last integral is computed via Eq. (3.17), and the final result is 

t~- 2 \ /,x r(2n-d/2) „ JJfel 2 , 

(*) = Ur(2n) Wd( 4 } ' ( } 

with F n d as in (1.15) or (1.16) (to obtain this, one also uses Eq. (1.13) for T). 
Inserting (4.13) into (4.11) we get the thesis. o 

Proof of Prop 2.2, item i). Lemmas 4.1 and 4.2 give immediately the bound 
(2.2) for K nd , with as in Eq. (2.3). 

We now pass to the boundary values of the function y nd for u = and u — > +oo. 
To determine =5^(0), use either Eq. (1.15) or Eq. (1.16), together with Eq. (3.7); 
the result agrees with Eq. (2.4). 

To determine lim u ^ +OQ S^ nd (u) we use Eq. (1.16), the limits 

' + l " V ' "■ " 1 forw^+cx), (4.14) 



1 + u J u + 1 

and Eq. (3.7); these relations imply 

_ 2" r(n+l/2)r(n-d/2) _ T(n-d/2) 
nA+ ' 7r d / 2 +V2 r(2n) 2^% rf / 2 r(n) ' 1 } 

where the last equality follows from (1.13). This gives the expression in (2.4) after 
using (1.10) with w = n — d/2. 

Of course, the continuity of y nd on [0, +oo) and the finiteness of its u — > +oo limit 
ensure that 5^ n d is bounded on its domain. o 

Proof of Prop 2.2, item ii). Step 1. The function y nd is increasing if d/2 < n ^ 
d/2 + 1/2. To prove this, we use Eq.s (2.3) (1.16) and the following remarks: 

a) the functions u G [0, +oo) — > (1 + 4w)/(l + u) G [1,4) and u G [0, +oo) -> 
u/ (1 + u) G [0, 1) are increasing; 

b) the function w G (— oo, 1) — > F(n,d/2 + 1/2 — n,n + 1/2; w) is increasing for 
d/2 < n < d/2 + 1/2, due to (3.13); in the limit case n = d/2 + 1/2, this function 
equals 1 everywhere (by (3.9), with m — 0). 

Of course, the fact that ,y nd is increasing implies sup[ 0j+oo ) S^ nd = y nd (+oo), and 
this fact, with Eq. (2.4), yields Eq. (2.5). 

Step 2. The asymptotics (2.6) of K+ d for n -> (d/2)+. This is evident from (2.5). o 

Now we must prove item iii) of the same proposition, concerning the n — > +oo 
behavior of K^ d ; a fairly long series of Lemmas will be established to this purpose. 
A main point in this argument is the integral representation, coming from Eq.s (2.3) 
(1.15) and (3.10), 

_ r(2n-d/2)r(n+l/2) 

" 2 ^/ 2 +V2r(n)r(2n) ' (416) 
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^( u ):=(l + 4 U )-/d. 7r =jr__. 

For future convenience, we write 

tf nd (u) = Kd(u) + <% nd (u) , (4.17) 



<,(„) := (1 + 4»)» / d. - ?r ^— ^ , (4.18) 



/•1/4 n-1 

^(m) := (l + 4 M ) n / rfs 



/ v / r^s(i + MS)" /2 ' 



4.3 Lemma. Define 
then, for fixed d and n — > +oo ; 



B nd := sup ^ nd (w) ; (4.19) 
«e[o,+oo) 



n 

Proof. We will estimate 38 nd (u) with different methods for u e [0, 2] and -u e 
(2, +oo), respectively. 

Let ^ u ^ 2; we reexpress the definition of S$ n d(u) as 

/•1/4 d/4-1 / N«-d/4 

*-<«> = a + *0"j( *7r^((TT^) ' (421) 

The function s — > s/(l + ws) 2 is increasing for ^ s < 1/u; but 1/u > 1/4, so the 
maximum of this function for ^ s ^ 1/4 is attained at s = 1/4. From here one 
gets 

/•1/4 d/4-1 / 1/4 \n-d/4 

+ "•^((I+W) = (422> 

ra-d/4 -1/4 „d/4-l 

On the other hand, the function w — > (1 + 4tt)/(2 + u/2) 2 is increasing for ^ u ^ 2, 
and equals 1 when w = 2; from here and from (1 + 4-u) d//4 ^ 9 d//4 one easily obtains 



/•1/4 s d/4-l 

sup @ nd {u) < C d for all n > d/2, C d := 9 d/4 / ds —== . (4.23) 
e[o,2l Jo V 1 - s 



ue[o,2] 
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We pass to bind £$ nd for u G (2, +00). Returning to Eq. (4.18), we write 



2 r 1/4 



,n-l 



2 (\ + Au) n r /4 



+ us) 2n - d / 2 



n-l 



(4.24) 



U n J (1 + g)2«-d/2 ' 

where the first inequality follows from — s < 2/\/3 for ^ s ^ 1/4, and the 

subsequent equality is obtained putting s = q/u. On the other hand, (1 + 4u)/u < 
9/2 for u > 2 and Jjj 1 ^ 4 < J +o ° on positive functions, so 



2 9 / >+ °° (/ n_1 

^,^ (tt)< ^ ( 2 )B y ^(tt^ 

2 f 9. n T{n-d/2)T(n) 



(4.25) 



v^V r(2n-d/2) 

(recall Eq. (1.14)). We now apply the duplication formula (1.13) with w = n — d/4; 
this gives 



sup &nd{u) < 
«e(2,+oo) 



Tin) 



r(n-d/4 + l/2) 



2 2+d/2 /9 Nn r(w-rf/2) 

V3 V r(n-d/4) 

Putting together Eq.s (4.23) (4.26) we get 

«^-( C '4 2W,, 0"^Hr ( »- i/4 + W 



(4.26) 



Tin) 



(4.27) 



for all n > d/2. As a final step, we recall that ([11], page 119) 



r (^ + fl ) =w a- b 



l + 0(- 

W ' J 



r(n) 



= 71 



-1/2 



r(w + 6) 

this implies, for n — > +00, 

T(n - d/2) 
r(n - d/4) r(n - d/4 + 1/2) 

and Eq.s (4.27) (4.29) yield the thesis (4.20). 
4.4 Lemma. For all n > d/2 one has 



sup £^ nd {u) < A nd , 
ue[o,+oo) 



for fixed a, 5 G R and u> — > +00 ; 



i + o(-) 

n 



(4.28) 



(4.29) 



o 



(4.30) 
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A 



nd 



->2n 



_ d/2 (1 - d/2n) 



i-d/2 pi 



(1 - d/4n) 2n - d ' 2 A /4 sVT^ (4 - s)«- d / 2 
For /ixerf d and n — > +oo ; z£ is 



-4nd = V^ 



3 d/2+l/2 4 



i + o(-) n 

n 



2 d / 2 ^ v 3' 

Proof. Step 1. The hound (4.30). The definition of &/ nd implies 



(4.31) 



sup s/ nd (u) < 
«e[o,+oo) 



'1/4 VI - S 



H nd (s) , 



(4.32) 



(1 + 4u) n 



H nd (s) := sup — , 2n _ d/2 ■ 

«e[o,+oo) (1 + su) 2n d ' 2 



For s E (1/4, 1), the function u E [0, +oo) -> (1 + 4w) n /(l + su) 2n ~ d/2 attains its 
maximum when u equals 



Und(s) ■-- 



8n+ (d-4n)s 
4(2n - d)s 



(4.33) 



Thus 



H nd (s) 



[1 + 4u) n 



(1 + su) 2n ~ d / 2 



u=u nd (s) 



d_^n-d/2 2 2n - d / 2 

(]_ d_y n -d/2 s n(4 _ g^n-d/2 



(4.34) 



and inserting this equation into (4.32) one gets the thesis (4.30) . 
Step 2. The asymptotics (4.31). We reexpress Eq. (4.30) for A nd as 



U nd ■ — 



A nd = 2 2n ~ d / 2 U nd f ds e(s)e 

Jl/4 

(1 " l) n ~ d/2 



•(ra-d/2)S(s) 



(4.35) 



I 4_\2n-d/2 ' 

An ' 



Q(s) := 



S\/l — s ' 



$(s) := log(4 - s) . 



In this representation we recognize a Laplace integral in the parameter n — d/2; it 
is $'(s) < for all s E (1/4, 1), $(1) = log 3, and the scheme of Eq.s (3.23)-(3.28) 
suggests to rephrase Eq. (4.35) as 



A nd = 2 2 - d / 2 U nd e~^ d ^L(n -±) = (^/ 2 (|)» U nd L(n - ±) , 



L(m) :- 



-3/4 



dt ${t)( 



-m<p(t) 



(4.36) 
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0(f) := 0(1 -t) = -^—) ' Vit) := $(1 ~ t] ~ = bg(1 + \ ] ■ 

The last integral has the standard Laplace form (3.18), and the framework of Prop. 
3.1 prescribes to analyze it introducing the function 

m := m = 3+« 3 

w ^(t) (l-t)Vt 

For t — > + , one has 

y?(t) = ^ + 0(t 2 ) , t = 3y2(t) + 0(y?(t) 2 ) , 

£(f) = ^ + O(Vt) = + O(v^)) • (4.38) 

Now, application of Prop. 3.1 to the last relation (4.38) gives 



'37T 1 

L(m) = + 0{— -pr) for m -> +oo . (4.39) 



On the other hand (taking the logarithm and expanding) 

U nd = 1 + O(-) forn -> +oo ; (4.40) 

n 

inserting Eq.s (4.39) (4.40) into (4.36), one easily derives the thesis (4.31) . o 



4.5 Lemma. For fixed d and n — > +oo ; it is 

.1 3W 4 

'2 } ~ Vn 2 d / 2 ^ ( 3 

(nofe f/ia£ f/ie n<//i£ /iand szrfes of this equation and (4-31) coincide). 

Proof. The definition (4.16) gives 



K d O = Gr 1 + °(-)} ( 4 - 41 ) 



1 Z" 1 S n_1 Z" 1 

-) = 3 n / rfs = 3™ / 

2 ; Jo v/r^(i + s /2)"/ 2 y 



rfs 9 (i (s)e- (n - d/4) * (s) , (4.42) 



s d/4-l 

Q d (s) := -== , $(s) := 2 log(l + s/2) -logs. 
V 1 - s 

We have again a Laplace integral, with parameter n — d/4; one finds &'(s) < for 
all s G (0,1), $(1) = 21og(3/2) and referring again to the scheme (3.23-3.28) we 
reexpress (4.42) as 

= re-WWL d (n -\) = (^ 2 (- 3 T L d (n - \) , (4.43) 
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L d {m) := [ dt $ d (t)e 
Jo 



rrup(t) 



(I _ tY^" 1 t 

Mt) := Q*(l-t) = , <p{t) := $(l-t)-$(l) = 21og(l--)-log(l-t) . 

Following again the scheme of Prop. 3.1, we introduce the function 

Mt) (3-t)(i-t)^ 
U >- v'it) (i + tyt ' 1 ' 

Let us keep d fixed. It turns out that Eq.s (4.38) are again satisfied with the present 
choice of ip and with £ = £ d . Therefore, Prop. 3.1 gives the asymptotics, analogous 
to (4.39), 



V3tt 1 

L d (m) = —= + O(^rrr) for m -> +oo ; (4.45) 
inserting Eq. (4.45) into (4.43) we obtain the thesis (4.41). o 

4.6 Lemma. For fixed d and n — > +oo ; one has 

_ Qd/2+1/2 A r 1-, 

sup = v^F -^-= (-)» 1 + O(-) (4.46) 

«e[o,+oo) z 7 v n <J L ^ J 

(again, the right hand side is as in Eq. (4.31)). 
Proof. We have 

^nd(-) < sup Kd(u) < A nd + 5 nd , (4.47) 

^ «G[0,+oo) 

(the upper bound follows from Eq.s (4.17) (4.19) and (4.30)). Both the above 
bounds on sup^ n d have asymptotics as in the right hand side of Eq. (4.31). For 
the lower bound, this is granted by Lemma 4.5. For the upper bound, this follows 
from Lemmas 4.4 for A nd and 4.3 for B nd : the latter is negligible with respect to 
the former, since 

B nd = A nd = A nd 0((j|r) = A nd O(-L) for any real a . (4.48) 



o 



Proof of Prop 2.2, item iii). Eq. (4.16) and the definition of K+ d in Eq. (2.2) 
give 



9> 1 / r(2n-d/2) r(n + l/2) / 1 

^nd{ 2 ) 2 , /27rd/4+1/4 y r(2n) r(n) y^nd{ 2 ), I4.4yj 
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1 Y{2n-d/2) T(n + l/2) 



«i = 2*v/«/ 4 y-r(2^ rw~ ylJ^^M ' (450) 

We know that ^^(1/2) and sup[ 0+oo ) ^ n( 2 have the same asymptotics up to 0(l/n), 
given by Lemmas 4.5 and 4.6; furthermore, Eq. (4.28) implies 



r(2ra-d/2) 1 r. _,l.i T(n+l/2) 



T(2n) (2n) d / 2 



l + 0(-) 
n 



l + 0(-) 

n 



(4.51) 



and inserting these results into Eq.s (4.49) (4.50) we obtain the thesis (2.7). o 

Computing the upper bounds K^ d . a) For d/2 < n ^ d/2 + 1/2, we have 
for the explicit expression (2.5); this was employed to compute the numerical 
values reported in Table 1 for these cases. 

b) In all the other cases, to compute K^ d one has to maximize the function y nd 
given by Eq. (2.3), containing the hypergeometric function F n d of Eq.s (1.15-1.16). 
For n — d/2 — 1/2 integer, 5^ n d has the elementary expression (1.17). 

c) Apart from simple exceptions, the maximization of S^ n d must be performed nu- 
merically. In all the cases analyzed with n > d/2 + 1/2, we have found numerical 
evidence (and sometimes an analytical proof) that y n d has a unique maximum point 
u = u n d > 1/2 in the interval (0, +oo), so that 

Kd = V^ndiUnd) . (4.52) 

d) Let us consider, for example, the case d = 2. For n = 2, Eq.s (2.3) (1.15) give 

one finds numerically that =5^2 attains its maximum at w 2 2 — 6.84. For n = 5/2, 
using (2.3) (1.17) one finds 

^ {l + Auf/ 2 Q + u (A KA . 

^ 2(M) = ^^(TT^ ; (454) 



the point of absolute maximum of this function is ^5/2,2 — 16/5 = 3.2, determined 
analytically by solving an algebraic equation of second degree. For larger, half- 
integer values of n, y n2 is again elementary, but the analytic determination of its 
maximum point involves algebraic equations of order increasing with n; thus, a 
numerical attack is necessary. 

Table 1 also considers, for d = 2, the values n — 4, 7, 16, 31, 61, 121. In all these 
cases, one finds numerically a unique maximum point u n 2 — 1.46, 0.915, 0.654, 
0.576, 0.538, 0.519. Note the approach of this point to the limit value u — 1/2 for 
large n, as expected from Eq. (2.7); due to this behavior, numerical maximization 
is simple even for very large values of n. 
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5 Proofs for the Bessel lower bounds on K n ^ 



Proof of Prop. 2.3, item i). Eq.s (2.11)-(2.12) are obvious; we must justify the 
expressions (2.13-2.14) of I^awIU, and (2.15-2.16) for ||^ nd || n . 
Step 1. Computation of \\g\ n d\\n- We have 

p i p (-\ _i_ \h\ 2 \ n 

hu\i = 1 dk (i + i*m w = ^ J Rd d k {l \^)l 2)2n = (5.D 
_ 2 ^ 2 r°° dPpd -i (!+p 2 ) n _ * d/2 r duu ^ g+^r 

"r(rf/2)A M 7 ^ (l + p 2 /A 2 ) 2 ™ T(d/2)\ d J {l + uY n ' 

In the last two passages we have used Eq. (3.1) for the integral of a radially sym- 
metric function, depending only on p := \k\, and then we have changed the variable 
to u = p 2 /\ 2 . 

For n arbitrary, the last integral in u is computed using the identity (3.11); this 
gives the thesis (2.13) (after using Eq. (1.10) with w = n — d/2). 
For n integer, in the integral over u we expand (l + \ 2 u) n with the binomial formula, 
and integrate term by term; this gives Eq. (2.14) after treating each term by (1.14). 
Step 2. Computation of \\gl nd \\ n - According to the definition (2.9), the function 
gxnd is obtained from the g n d of Eq.(1.20) rescaling by A. From here, and from Eq. 
(4.13) for Tg 2 nd we infer 

W )(h) 1 nin-d/2) \k\ 2 

[J-g Xnd ) {k) - 2d/2xd v{2n) *nd{ 4X2 ) , (5.2) 

with F nd as in Eq.s (1.15) or (1.16); thus, 

\\9l d \\l= [ dk (1 + \k\ 2 r \Tg 2 Xnd (k)\ 2 = (5.3) 

_ ^ 2 T 2 (2n - d/2) [+» p 2 
~ 2^T(d/2)T*(2n)\™ J dp P (1 + P ] F «*Up " 

Now, introducing the scaled variable u := p 2 /(4 A 2 ) we readily obtain the expression 
(2.15) for \\glJ\ n . 

Finally, let us consider the case n — d/2 — 1/2 integer and show that Eq. (2.15) 
becomes Eq. (2.16). In fact, in this case the function F nd has the elementary 
expression (1.17); when this is substituted into the integral over u of Eq. (2.15), we 
get 

/ duu d/2 -\l + 4X 2 u) n F 2 d {u) = (5.4) 
Jo 

n ~ d ^ 2 ( W ), (d/2 + 1/2-71), (n) m (d/2 + l/2-n) m 
e ^ (n+l/2) e £\ (n + l/2) m m! 



27 



y (i + u ) 2n+e+m 

each of the above integrals can be computed via Eq. (3.11), and the conclusion is 
the thesis (2.16). o 

To prove the second item in Prop. 2.3 we need an elementary bound for the 
hypergeometric-like function F n d, to be substituted in Eq. (2.15) for ||<7> n J|n; this 
will require some Lemmas. 

5.1 Lemma. Assume 

f G C([0, 1], R) n C 2 ([0, 1), R), ReC([0, l],R)f]C 1 ([0,l),R), e>0; (5.5) 

f(w) = (1 - wY^Riw) , R'(w) > for we [0, 1), (5.6) 
and consider the C 2 function 

Then: 

/( i)-/H J(r< ^ r 



(1-wf e 
d /(1)-/H 



(1 — w) ( 

The previous facts imply 



> for we [0, 1) . (5.9) 



/(1)-/(0)<M^M<M /or «, e (0,1). (5.10) 
(1 — w) £ e 

Proof. By the generalized Lagrange theorem, it is 

F(l) - F(w) F'(t w ) 



G(l) - G(w) G'(t w ) 



for some t w G (w, 1) , (5-H) 



for all F,G G C([0, 1],R) n C\(0, 1),R) with G" never vanishing, and for all w G 
[0, 1). We apply this statement with 

f:=F, G(w) := -(1 - wY , (5.12) 

taking into account Eq. (5.6); this gives 

/( (i-i) w = f ° r w e [0, i}> with tw g {w > 1)5 (s - i3) 
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and in the limit w — > 1 we obtain Eq. (5.8). 
In order to prove (5.9), we observe that 

f{ ]] " f [ W) = 7r -i TTr (/(l) - /H) - ^ for W G [0,1). (5.14) 
rfw (l-w) e (1 - w v ' y 11 l-w 17 v y 

On the other hand (intending as an improper Riemann integral) 

e (/(l) - /H) = ef dt f\t) = e f dt(l- tf-'Rit) = 

= (1 - w) e R(w) + I dt(l-t) e R'(t) . (5.15) 

J w 

the last equality following from integration by parts. Inserting (5.15) into (5.14) we 
obtain 

d /(i)-/H i 



dw (1-wY (l-w) e+1 



j dt(l-t) £ R'(t) , (5.16) 

J w 



and the positivity of R! gives the thesis (5.9). 

Finally the function w G (0,1) i— > (/(l) — /(u>))/(l — u>) £ is increasing, so it is strictly 
bounded from below and above by its limits for w — > + and w — > l - ; this yields 

Eq. (5.10). o 

5.2 Lemma. Let 

0<a,b<+oo; a + b < c < a + b + 1 ; w e (0,1) . (5.17) 

r>/ 7 \ P(a,b,c) - F(a,b,c;w) , , /rio > 

< P(a, 6, c) - 1< -5= — 7 _ ^ c _ a _ h < Q(a, b, c) (5.18) 

where 

P(a,6,c): = F(a,6,c;l) = r ; c ^ a)r(c _ 6 j , (5.19) 



(c - a - 6)r(a)r(6) 
Proof. We apply the previous Lemma with 

f :=F(a,b,c;.) , e:=c-a-b. (5.20) 

In this case, the differentiation formula (3.12) and the subsequent application of the 
Kummer transformation (3.15) give 

db 

f(w) = (1 - wY^Riw) , R(w):=—F(c-a,c-b,c + l;w). (5.21) 

c 
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On the other hand, the hypergeometric function w i— > F(c — a, c — b, c + 1; w) has 
positive derivative, due to (3.13) and to the assumptions (5.17) for a, b, c; the same 
assumptions ensure this function to be continuous also at w — 1, where its value is 
determined by Eq. (3.7). Thus all conditions of the previous Lemma are fulfilled by 
f,e,R, and Eq. (5.10) gives 

t?( h i\ t?( h n\ / F(a,b,c; 1) - F{a,b,c;w) 
F(a, b, c; 1) - F(a, b, c; 0) < (1 _ w y- a - h < 

<-. — -F(c-a,c-b,c + l;l) . (5.22) 

c(c — a — o) 

But 

F(a,b,c;l) = P(a,b,c), F(a, b, c; 1) - F(a, b, c; 0) = P(a, b, c) - 1 > ; (5.23) 

the last inequality holds because F(a,b,c;-) is increasing (see again Eq. (3.13)). 
Finally, the equality 

F(c-a,c-b,c+l;l) = Q(a,b,c) (5.24) 



c(c — a — b) 

is easily inferred from Eq. (3.7), using the identity (1.10) with w = a and w — c. 
Eq.s (5.22)-(5.24) yield the thesis. o 

Remark. The idea of employing (3.15) in the above proof has been suggested by 
[12], where the usefulness of this transformation has been pointed out in relation to 
similar inequalities for F. 

5.3 Lemma. Leta,b,c, P(a,b,c),Q(a,b,c) be as in Lemma 5.2, and 

q(a, b, c) := { Q ^ & ' C) ^ ^ ^ ^ ^ Q ^ ^ ^' (5.25) 

| P(o,ii,c)-l if P(a,b,c) < Q(a,b,c). 

Then 

F(a, b, c; wf > P(a, b, cf - 2P(a, b, c) Q(a, b, c)(l - u>) c ~ a ~ 6 + (5.26) 

+q(a, b, c) 2 (1 - w) 2{c - a - b) for w G (0, 1). 

Proof. Step 1. The case P(a, b, c) ^ Q(a, b, c). For any w G (0, 1), the upper bound 
in Eq. (5.18) implies 

F(a, b, c; w) > P(a, b, c) - Q(a, b, c)(l - wf- a - b . (5.27) 
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The right hand side in the above equation is positive, so we infer 

F(a, b, c; wf > (P(a, b, c) - Q(a, b, c)(l - w) c - a ~ b ) 2 ; (5.28) 

expanding the right hand side we get the thesis (5.26), since in this case Q(a, b, c) = 
q(a, b, c). 

Step 2. The case P(a, 6, c) < Q(a, b, c). We write 

F(a, b, c; wf = [P(a, b, c) - (P(a, b, c) - F(a, b, c; w)) ] 2 = (5.29) 

= P(a, b, c ) 2 + (P(a, b, c) - P(a, b, c; w)) 2 - 2P(a, b, c)(P(a, b, c) - P(a, b, c; w)) . 

We insert here the bounds on P(a, b, c) — P(a, b, c; w) coming from Eq. (5.18); this 
gives 

F(a, b, c; w) 2 > P(a, b, c) 2 + (P(a, 6, c) - 1) 2 (1 - w) 2(c ~ a " fe) + (5.30) 

-2P(a, b, c)Q(a, b, c)(l - w) c ~ a ' b , 

and we have the thesis (5.26) since in this case q(a, b, c) = P(a, b, c) — 1. o 

Proof of Prop. 2.3, item ii). Throughout the proof, d/2 < n ^ d/2 + 1/2. 
Step 1. For w G (0, 1) one has 

d 1 1 



F(ra,- + --n,ra+-;u;J ^ (5.31) 

^(n-d/2)=» (n-d/2)2 U Wj + („ - d/2)^ 1 Wj ' 

where P nd) Q nd and are as in (2.21). For n < d/2 + 1/2, this follows from 
application of Lemma 5.3 with a — n, b = d/2 + 1/2 — n, c = n + l/2; comparing 
the coefficients in this Lemma with Eq. (2.21) we see that 

p K b i c ) = Pnd A ,~ » <2( a ' b > c ) = ^ n j/n > 0( a » & > c ) = /0 • ( 5 - 32 ) 

n — a/2 ra — a/2 n — a/2 

Let us pass to the limit case n = d/2 + 1/2; then, (5.31) holds as an equality 
because P nd = 1/2, Q nd = 0, q nd = 0, F (n, d/2 + 1/2 - n, n + 1/2; w) = F(d/2 + 
1/2, 0, d/2 + 1; w) = 1 (by (3.9), with m = 0). 

Step 2. Proof of Eq. (2.17): \\g 2 Xnd \\ 2 n > S?„d(A), with <g nd (\) as in Eq. (2.20). We 
start from the expression (2.15) of ||g 2 n(i || 2 ; the function F nd therein is expressed as 
in (1.16), and its square is bounded via the result of Step 1 (with w — u/(l + u)). 
This gives 

Ik 2 II 2 > - d/2 ^n-d/2) 

\\9xnd\\n * {n _ d/2 yr(d/2)ri(2n)\i x 
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X 



r duu d '^ {i+Ax2u)n (p* 2 PndQnd i 



Qnd 



+ u y n -dj 

The above integral can be written as the sum of three integrals of the form (3.11); 
after computing each of them by (3.11), we apply (1.10) with w = n — d/2, 2n — d 
and 3n — 3d/2, respectively. The final result is the minorant for ||g 2 nrf || 2 as in Eq. 
(2.20). 

Step 3. The n -> (d/2)+ limit of Jt n B d B {\). Let d and A G (0, +oo) be fixed. We 
start computing the limiting behavior of ^^(A). For n — > (d/2) + , the coefficients 
P n d-> Qnd and q n d therein have the same behavior up to 0(n- d/2): 



P O a r ( rf / 2 + 1 /2) 

r ndi Wndi Hnd 



l + 0(n-d/2) . (5.34) 



v^r(rf/2) 

In the same limit, the three hypergeometric functions also have equal behavior: 
F(-n,^,n;l-4A 2 ), F(-n, ^ 2n - ^ 1 - 4A 2 ), F(-n, ^ 3n - d; 1 - 4A 2 ) = 

= F{-± d - d -;l- 4A 2 ) + 0(n -±) = 2 d \ d + 0(n - ±) , (5.35) 

where the last equality follows from (3.8). Inserting Eq.s (5.34) (5.35) into (2.20), 
we find 

2 V/ 2 - 1 r(rf/2 + 1/2) 2 1 + Q(n - d/2) 
nd( j ~ 3T 2 (d)r(d/2) (n-d/2)3 

^ l + Q(n-d/2) . 
3 2 d - 2 r(d/2) 3 (n-d/2) 3 ' 1 } 

the second equality in (5.36) follows from the first one applying the duplication 
formula (1.13) with w = d/2. 

Let us pass to the n — > (d/2) + behavior of || (?A«d||n; from (2.13) and (3.8), we infer 
II^awL - r{d/2)Xd n-^ 2' 2' 1 " A } n -d/2 ~ (5 - 37) 



7T 



d / 2 l + 0(n-d/2) 



r(d/2) n-d/2 



Since JC B (A) = v/CKAl/lkAndll 2 , from (5.36) and (5.37) we obtain 

^(A) - ; - \ , 1 + °^~ d ^ ; (5.38) 

V ; v / 3 2 d / 2 - 1 7r d / 4 v ^72) - d / 2 

comparing this with the definition (2.6) of M d , we get the thesis (2.22). o 
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Computing the Bessel lower bounds, a) For all n > d/2, the lower bound 
J^J(A) is the ratio of llgfndlU an d ||<7And||n- The norm of g\ nd has the analytic 
expression (2.13) in terms of a hypergeometric function, that becomes the elementary 
formula (2.14) for n integer. 

The norm of g\ nd has the integral representation (2.15), involving the hypergeometric- 
like function F nd of Eq.s (1.15) (1.16). For n — d/2 — 1/2 integer, this norm has the 
explicit expression (2.16) in terms of hypergeometric functions. For n — d/2 — 1/2 
noninteger, the integral in (2.15) must be computed numerically. As anticipated, 
this is a difficult task for n very close to d/2, due to the slow convergence of the 
integral: the integrand behaves like l/u 1+ ^ n ~ d ^ for u — > +oo (as made evident by 
Eq. (1.16) for F nd ), and we are interested in situations where n — d/2 = 10~ 4 . In 
these cases it is convenient to compute, in place of Hs'Andlln) the minorant Sf n <i(A) of 
Eq. (2.20), and from this the lower bound X^ B (X) of Eq. (2.19), both of them 
having analytic expressions in terms of hypergeometric functions. 

b) Assuming we are able to compute J^(A) or J^ S (A), for each A we have a lower 
bound for K nd ; the next step is maximization with respect to A, to get K^ d or K^f . 
In general, this is done numerically (using some package for automatic maximization 
or for plotting these functions of A, so as to read the maximum from the graph). 

c) Let us consider, for example, the case d = 2 and the values of n reported in Table 
1. For n = 3/2, we have the elementary expression 

F(w) := F(-3/2, 1, 3/2; w) = + |(1 - w) 2 ^(w) , (5.39) 

arctanh( v /w)/ v /ii; if < w < 1 , 
&(w) := { 1 if w = , 



arctan(A/— w) / ^/—w if w < . 

The function ^^ 2 2 stains its maximum at A ~ 1.38. ^^ 2 2 lii a ^ so elementary, 
with its maximum at A ~ 1.36. For n — 5/4, 2, 4, 7, 16, 31, 61 the integral in ||^ n2 IU 
can be computed numerically; from the graph of JXf^ we have found this function 
to get its maximum at A ~ 1.40, 1.36, 1.39, 1.45, 1.53, 1.57, 1.58, respectively. 
For n — 1 + 10~ 4 , 1 + 10~ 2 , 1 + 10" 1 the numerical computation of J^f and K^ 2 
is difficult, so we have turned the attention to the simpler bound K®®\ from the 
analytic expressions of ^^f B (A) and numerical optimization, we have found the 
maximum of this function to be attained at A ~ 1.42 in each one of the three cases. 
For all the cases in the table from n = 5/4ton = 61, the previously mentioned 
Bessel bounds have been compared with the Fourier lower bounds K^ 2 or K^ 2 F of 
Prop. 2.4 (for the computation of these Fourier bounds, see the remarks at the end 
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of the following section). In this way, we have found that the Fourier lower bounds 
are below the Bessel bounds up to n = 5/2, while the contrary happens forn > 3/2 
(for example: #jL 2 < 0.610 iriw 4)2 and Kg 2 < 0.411 K^ F 2 ). Extrapolating, the 
Bessel bound K^ 2 is likely to be smaller than the Fourier bounds for the large value 
n = 121. Since the numerical computation of J^fi j2 and Kf 212 is difficult, in the 
construction of Table 1 we have chosen directly for K^ 212 a Fourier bound. 



6 Proofs for the Fourier lower bounds on K na >. 

We refer to the trial functions f pad of Eq. (2.23). Our aim is to prove all statements 
contained in Prop. 2.4; we will proceed in several steps. 

6.1 Lemma. For all p,a > and n > d/2, \\f pa d\\n is given by Eq. (2.26). 
Proof. The Fourier transform of f pad is elementary, and given by 

(Ffv*d)(k) = ^e-£l fe -H 2 , v := (1,0,. ..,0) ; (6.1) 

thus 

\\f P «Al = - d [ dk(l + \k\ 2 Te-^\ 2 = (6.2) 

= \f dk(l + \k\re-^ + ^ k . 
To go on, let us first consider the case d — 1. Eq. (6.2) gives 

\\f P ,i\\l = - [ dk(l + kre'^^ k = (6.3) 
c Jr. 



dp(l + p"re-^(e^ + e-^) 



(in the last passage, we have used the variable p = \k\); this gives Eq. (2.26) for 
d = 1, since ([14], page 80) 

e s + e~ s = V2^ I-i/ 2 (s) Vs G (0, +oo) . (6.4) 

Now, let us pass to the case d ^ 2. Eq. (6.2) contains an integral of the form (3.2), 
where the integration variable is now k and x(|&|> V*k) = (l + |£;| 2 ) n e~^' fe ' 2+p2 ) + ~^' ?,fe ; 
therefore, 

9 d/2-1/2 r+oo 2 , 2 r n i 

\\f P « d \\l = r(d n d / d P P d -\1 + P 2 ) n e-^ / d9 sinO^e^ . (6.5) 
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On the other hand ([14], page 79) 

£ d6 sm6 2u e scosd = v^r(z/+ i) 0^ I u {s) ; (6.6) 
inserting this result into the previous equation, we obtain the thesis (2.26). o 

6.2 Lemma. For allp,a > and integer n > d/2, \\f pa d\\n is given by Eq. (2.27). 

Proof. We return to the first equation (6.2), and expand (1 + \k\ 2 ) n by the binomial 
formula; this gives 

Now, we write the integration variable as k = (h, q), (h G R, q G R^ 1 ); so, 



ll/Jn = iE( ? / _ dhdq(h> + \q\ye-^e-^= (6.8) 



2 ' 9 | 2 



a {Zg V I I JnxR d - 



1 x / n 



J e ct 



where, in the last passage, we have used again the binomial formula to expand 
{h 2 + \q\ 2 ) e - On the other hand, 



L 



dhh 2j e~ (Jl ^ = [ dh(h + p) 2j e~^ = (6.9) 



Uiy "I.-"--* -w.) 



2j \ {2g-^Wl 1/2+g 



-p J "a 
29 y 



The last passage above depends on the evaluation of the integrals with h m : these 
vanish for m odd, while in the even case m = 2g we have f + °° dh h 2g e~ h2 ^ a = 
a 9+1 / 2 r(g +1/2) = a 9+1 l 2 2~ 9 (2g — 1)!!^- Concerning the integrals over q, due to 
Eq. (3.1) we have 

P 2 O d/2- 1/2 P+OO 2 

J d i dq\qf- 2 >e-^ = J* ln 1 ^ / £+-M-V e -*r - 



V(d/2 - 1/2) 
!/2_ 

T(d/2 - 1/2) 



d/2-1/2 d/2-l/2+l-j r(d/2 — 1/2 + t — j) 
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= -K d l 2 - l l 2 a d ^ 2+l ~^d/2 - l/2) t .j . (6.10) 

Inserting Eq.s (6.9) (6.10) into (6.8), we finally get the thesis (2.27). o 

Proof of Prop. 2.4, item i). This is given by the two previous Lemmas. o 

We pass to item ii) of the same proposition, whose proof is more lengthy. The initial 
step concerns the expression of ||/ po -d||n when p is arbitrary and a = c/n (c > 0); in 
this case, the already proved Eq. (2.26) becomes 

9 d/2 (1/2+1 r+oo 2 2 9 

II WE = %*^r I ^/ /2 (l + pYe-"^W^p). (6.11) 

We will analyze this formula in the limit n — > +oo. In the first Lemma, p and c will 
be arbitrary; in the subsequent ones, based on the theory of Laplace integrals, we 
will consider a specific choice, ultimately yielding Eq. (2.29). 

6.3 Lemma. Fix p > 0, c > and d; for n — > +oo ; it is 

2 ^d/2-1/2^/2+1/2 
\\fp,c/n,d\\n = c d/ 2 +l/2 p d/2-l/2 [ X pc,d/2-l/2 (™) + (6.12) 

+o( Xpc4/2 - 3/2{n) )+o((i+ P r) , 

/•+OO _ 2 

X pca (n):= dpp a (l +p 2 ) n e- nJL ^ L - for all a E R . (6.13) 

Proof. We start from the relations 

I d/2 ^{s) = -g= h d {s) = (l + for all 8 G (0, +oo), 

6 d ,/i d eL°°((0,+oo),R) , (6.14) 

reflecting the asymptotic behavior of the Bessel functions I v {s) for s — > + and 
s — > +oo (see [14]). 

To go on, in Eq. (6.11) we write J + °° = f p + °° + JJf; in these two integrals, we 
substitute the representations (6.14) of Id/2-1 involving, respectively, b d and h d . 
This gives 

^d/2-1/2^/2+1/2 

ll/p.c/n.dll^ = ^2+1/2^/2-1/2 pW/2-i/ 2 (™) + r pcd (n) + Z pcd (n)\ , (6.15) 
where the X term is defined following Eq. (6.13), and 

Y pcd (n) ■= ^- [ +X ' dpp d l 2 -*l 2 {l + p 2 Ye-^b d {^p) , (6.16) 
2pn L c 
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J pcd 



n) 



c h d { p) 



We estimate these two integrals. Let B d : = sup( 0+oo ) \b d \, H d : = sup( 0+oo ) \h d \; then 

B d c 



Ypcd{n)\ ^ 



2pn 



-X, 



pc,d/2-3/2 



n 



(6.17) 



dpp d i 2 -^ = H d {i+ P 2 r^ 



4/2+1/2 
d/2 + 1/2 ' 



whence 



>W(n) =0( Xpc '" /2 ' 3/2(n) ) , Z^n) = 0((1 +p 2 )") for n -> +oo . (6.18) 



Substituting Eq. (6.18) into (6.15) we obtain the thesis (6.12). 
To go on, we observe that Eq. (6.13) can be rephrased as 

f+OO 



O 



^"pca(^) 



/ dpfe**^) , $ pc(p ) ;= i£ — ^- - log(l + p 2 ) 

Jp c 



(6.19) 



In the sequel, we apply the Laplace analysis to the integral (6.19). We will consider 
the special choice 

p:=^=, c=- A (6.20) 

and its double (2p, 2c): this makes easy to compute the minimum point of $ pc 
and <3>2 P ,2c- We repeat here the remark made in Sect. 2, after stating Prop. 2.4: 
different choices of (p, c) complicate the computations, with no sensible increase in 
the dominant term of the Fourier bound J(f F (p, c/n). (This conclusion is the result 
of a tedious analysis, that is not worthy to be reported here). 

6.4 Lemma. Let p,c be as in (6.20). For fixed a G R and n — > +oo ; it is 



X (n) = h^ ( 3 /2) 



2^/2+1/2 e n/6^ 
on 



n 



3 n/3 



1 + (Di- 



rt 



Proof. Step 1. Proof of Eq. (6.21). We put for brevity 

X a (n) := X pca (n) , $ := $ pc . 

Explicitly 

$(P) = ^(P-^ 1 ) 2 -I0g(l+P 2 ): 



(6.21) 
(6.22) 

(6.23) 
(6.24) 
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it is easily checked that 



\ 2. 1 M + V2p + 2 < < 1 

$ (p) = -(P ;=) — £ = for p = — , 



(6.25) 



Now, following the scheme of (3.30) we reexpress the integral under examination as 



XJn) = e 



(3/2) r 



■ nHl/V ~ 2) [L-(n) + L»] = [L-(n) + L+(n)] , 

•1/(2a/2) 



(6.26) 



/•1/(2V2) _ /-+oo 

L-(n) := / (ft ^(t)e-^W , L+(n) := / dt 

Jo Jo 



-ntp+(t) 



m)--=(^Ttr, (6.27) 

^(t) := T t) - H^) = T—t + J* 2 - log(l T — t + f) . 

The above two integrals have the standard Laplace form discussed in Prop. 3.1. 
Following the usual scheme, we fix the attention on the functions 



cTM .= = 1 3 T 2y^t + 2t 2 1 

For t — > + , one easily checks that 

^M = y**=F^f + 0(f), 



(6.28) 



(6.29) 



2 a/2 



1 



9 1 9a 3 
20t T 71 To ~ lCT 

3 ,9a 1 



+ <?(*) = 
+ 0( V / ^W) • 



(6.30) 



2 a/2+l/2 

We can now apply Prop. 3.1 to both integrals LJ(ra); this gives 



1 



2«/2+l/2 



3^ ,9a 1,1 
=F 



2^^ 10 5 nJ 



+ 0( 



n 



3/2- 



for n — > +oo 



(6.31) 
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and substituting these expansions into Eq. (6.26) we get the thesis (6.21). (Note 
the mutual cancellation of the terms =f(9«/10 — l/5)(l/n), in agreement with the 
remark concluding Sect. 3). 
Step 2. Proof of Eq. (6.22). In this case, we put 



One has 



X a (n) := X 2p:2c>a (n) , $ := $ 2 p, 2c • 



$(p) = ^(p-^ f ) 2 -log(i + P 2 ); 

^)^(p- ^ 2 ;f; +1 |o fo rp|^, 

$(>/2) = ^-log3. 



(6.32) 
(6.33) 



We can write 



nn 

X a (n) = e- n ^[L-(n) + L»] = [L"(n) + L+(n)] 



(6.34) 



/•1/V2 /-+00 

L-(n) := / ^(Oe-^W , L+(n) := / <ft 

Jo Jo 

0j(f):=(V2=F*) a , (6-35) 
^(t) := $(^2 =F £) — HV2) = T 2 -^t + I* 2 - log(l T^t+ \f) . 
We introduce the functions 

&?(t) 3 3^2y/2t + t 2 



-n<p+(t) . 



<p*'(t) 2t 7T$V2t + 2t 2 
For t — * + , comparing the expansions of ip T , £J in powers of t we get 



(6.36) 



£(t) = 2"/ 2 



1 ,9a 



2 V / 7 V / ^ ? W y/2 K U 49 y J 
Applying Prop. 3.1 to L^(n) we obtain 



(6.37) 



n 



2 a/2 



3v^F 1 9a 2 1 



rr 



3/2 < 



for n — > +00 , (6.38) 



and substituting these expansions into Eq. (6.34) we get the thesis (6.22). 



o 
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6.5 Lemma. Let p,c be as in (6.20). For fixed d and n — > +oo ; it is 



||/p,c/n,d||n 
||/2p,2c/n,d||n 



2 3d / 2 7r d / 2 (3/2) n nd/2 



3(4/2-1/2^/5 e™/ 6 



i + o(-) 

n 



3^2-1/2^/7 e ™/3 



d/2 



n 



(6.39) 
(6.40) 



Proof. To prove Eq. (6.39), we note that (6.21) implies 



_ 3y^75 (3/2)" 
A pCjC(/ 2-i/2lnj - 2d/4+1/4 e „ /6 ^- 



l + 0( 



(6.41) 



^cd/2-3/2(n) _ (3/2)" 1 
n e n ^y/n n' ' 

We insert these results into Eq. (6.12) for ||/ P)C / n) d||^, taking into account that the 
present choices of p, c imply 

3^/2+1/2 



c d/2+l/2 p d/2~l/2 



27rf/4+l/4 ' 



1 +P 2 



9 _ 3/2 
8 ~ 



9, 0.8 < 9 < 0.9 ; 



(6.42) 



e n/e " e n ^y/n ~ y n' 
The proof of Eq. (6.40) is very similar, depending on Eq.s (6.22) (6.12) . o 

Proof of Prop. 2.4, item ii). This item concerns the n — > +oo limit for the 

special Fourier lower bound K FF \ comparing the definition (2.28) of this bound 
with the notations of this section, we see that 



K 



FF 

nd 



ll/2p,2c/ 



n,d\\n 



\\fp,c/n,d\\n 

From Eq.s (6.39) (6.40) we infer, for n — > +oo, 



(p, c) as in (6.20) . 



(6.43) 



K FF = 

^nd 2 d 7T d / 4 



n 



d/A 



(6.44) 



l + 0(- 



n 



(5/3)V2 (2/> /3)n 



1 



ra J 



71/4 a n d/4 

In the last passage we have used the definition (2.7) of Td, our result is just the 
thesis (2.29). o 
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Computing the Fourier lower bounds, a) For any n and d, the function 
(p,cr) — > J^(p, a) in Eq. (2.24) is determined by the function (p,a) — > ||/ po -d||n- 
For n noninteger and given (p, a), this can be computed via Eq. (2.26), evaluating 
numerically the integral therein; for n integer, we have the elementary expression 



The bound K^ d is obtained maximizing J£T^(p, a) with respect to (p, a) e (0, +oo) 2 ; 
in typical situations this must be done numerically, even for integer n (in any case, 
the maximization problem is not dramatic because J£^(p, a) is a lower bound for 
all choices of (p, a), even not close to the maximizing pair). 

For very large values of n, instead of maximizing J^(p,a) one can evaluate it at 
(p,a) = (l/(2V^),3/(4n)); this yields the the bound KfJ of Eq. (2.28), that we 
know to be effective in this limit. 

b) Let us consider, for example, the case d = 2 and the values of n in Table 1. For the 
integer values n = 2,4, 7, 16, 31 we have determined the analytic expression of JtJ^ 
using Eq. (2.27), and then maximized this function numerically; the maxima occur, 
respectively, at (p, a) ~ (0.511,1.05), (0.417,0.309), (0.371,0.148), (0.331,0.0582), 
(0.316,0.0290). For the large values n = 61, 121, we have used directly the lower 
bound = ^n2(l/(2\/2), 3/(4n)); since n is integer, in principle this could be 
obtained again from Eq. (2.27), but in these two cases it is more convenient to 
compute it numerically, starting from the integral representation (2.26) of H/p^Hn 
(note that this contains the non elementary function I ). 

For n = 5/4,3/2,5/2, (p, a) has been computed numerically for many sample 
values of (p, a), starting again from (2.26); in this case, approximate maximization 
has been performed choosing the best value in the sample. The maxima are attained 
at p ~ 0.354 in the three cases, and a ~ 5.22, 2.41, 0.696, respectively. 
The numerical computation of J£^f and K% d is difficult for the small values n = 
1 + 10~ 4 , 1 + 10~ 3 and 1 + 10 _1 . On the other hand, for the reasons already 
explained at the end of the previous section the Fourier bounds should be below 
the Bessel bounds for these extreme values of n; therefore to construct Table 1 in 
these cases we have given up computing K^ 2 , and we have chosen directly for K~ 2 
a Bessel lower bound. 



A Appendix. The integral I^ u {h). 

This integral is defined by Eq. (3.16); we want to prove Eq. (3.17). We start from 
the identity ([14], page 440) 



(2.27). 




+oo 



(A.l) 
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and insert it into (3.16); this gives 

"+oo p+oo 



I^(h)=2 dt dr r^^ 1 J^(hr)K u (2r cosht) . (A.2) 
Jo Jo 

On the other hand ([14], page 410) 

' drr^+%(hr)K„(2rcosht) = ^l^t * ^ 

( h 2 \ r(u + v + i)h» ( h 2 x-"-"- 1 

xF[u + u + l,u + l,u + l; 5- = — — 9 , , 9 1 H 5- 

V 4cogh 2 t ; 2^+ 2 cosh 2/i+iy+2 1 V 4 cosh 2 

where the last passage depends on (3.8). Returning to Eq. (A.2) we obtain 

T r(/i + z/ + i)^ /-+ 00 , 1 / h 2 



(A.4) 



2»+ l in "cosh 2 ^ +2 tV 4 cosh 2 1 



2 M+2 



^ rf S S ^ 2 (l - S )-V2(l + ^S)-^" 1 , 



the last passage following with the change of variable s = 1/ cosh 2 1. Now, compar- 
ison with (3.10) gives the thesis (3.17). 



B Appendix. Proof of Prop. 3.1 on Laplace inte- 
grals. 

We recall the notations and assumptions (3.18-3.21), and point out some conse- 
quences of our hypotheses. 

First of all, by the monotonicity of tp, ip(b) := lhm,^- (p(t) exists in (0, +00], and <p 
is a C 1 diffeomorphism between (0,6) and (0, <£>(&)). 

Moreover, by Eq. (3.21), there are a constant e G (0,6) and a bounded function 
P G C((0,e),R) such that 

e-i 

£(t) = J2 P > ^T 1 ' 1 + MM*)*'- 1 f or all t G (0, e) . (B.l) 

i=0 

Putting the attention to Eq. (3.18) and dividing integration in two parts, we get 

L{n) = M(n) + N(n) , (B.2) 

M(n) := dt $(t) e- n ^ , N(n) := dt &(t) e~ nip(t) . 
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Let us estimate M(n). Introducing the new variable s = tp(t) and then using (B.l) 
we obtain 

M{n) = /" ( ) ds t{<p-\s)) e' ns = V P t M t (n) + 5M e (n) , (B.3) 



Mi(n) : = / ds s Q ^V ns , 5M e (n) := / ds s * -1 P{y-\s))e 
The above integrals are related to the incomplete Gamma function 

7(a,u):=/ dvv a - l e~ v = Y(a)- dvv a ~ l e~ v = 

JO Ju 

= T(a) + O {u a - 1 e~ u ) for u -> +00 (a > 0) 

(concerning the asymptotics of 7 for u — > +00, see [11]). As for M i: with a variable 
change s = v/n we get 



(B.4) 



7 («„^))^ + 



-n<p(e) 



n 



for n — > +00 ; (B.5) 



furthermore, 



\5M e (n)\ < (sup 



(sup 

(0,e) 



^ds S —e- = (sup'- ^^)) 

(0,6) 



(B.6) 



+ 



,—nip(e) 



n 



O 



1 



for n — > +00 



-(n— m)ip(t) p -mip(t). 



To estimate N(n), we fix ni > no and write jV(ra) = J^dt e 
for all n G [ni, +00), this implies 

|JV(n)| < e - (n - ni) ^ (e) y e"" 1 ^ = O (e- nv(e) ) for n -> +00 (B.7) 



(recall that < </?(e) < y?(t) for i e [e, 6)). 

From Eq.s (B.2), (B.5), (B.6) and (B.7) we get the thesis (3.22). 



o 



Acknowledgments. We acknowledge D. Bambusi and S. Paveri Fontana for useful 
bibliographical indications. This work has been partially supported by the GNFM of 
Istituto Nazionale di Alta Matematica and by MIUR, Research Project Cofin/2004 
" Metodi geometrici nella teoria delle onde non lineari e applicazioni" . 



43 



References 

[1] M. Abramowitz, LA. Stegun, Handbook of mathematical functions, Dover, New 
York (1992). 

[2] R.A. Adams, Soholev spaces, Academic Press, Boston (1978). 

[3] N. Aronszajn, K.T. Smith, Theory of Bessel potentials. I, Ann. Inst. Fourier 
11, 385-475 (1961). 

[4] S. Bochner, K. Chandrasekharan, Fourier transforms, Princeton Univ. Press, 
Princeton (1949). 

[5] Y.L. Luke, The special functions and their approximations, Academic Press, 
New York (1969). 

[6] V.G. Mazjia, Soholev spaces, Springer, Berlin (1985). 

[7] C. Morosi, L. Pizzocchero, On the constants for some Soholev imbeddings, J. 
Inequal. Appl. 6, 665-679 (2001). 

[8] C. Morosi, L. Pizzocchero, On the constants in some inequalities for the Soholev 
norms and pointwise product, J. Inequal. Appl. 7, 421-452 (2002). 

[9] C. Morosi, L. Pizzocchero, Quantitative functional calculus in Soholev spaces, 
Journal of Function Spaces and Applications (JFSA) 2, 279-321 (2004). 

[10] C. Morosi, L. Pizzocchero, On approximate solutions of semilinear evolution 
equations, Rev. Math. Phys. 16, 383-420 (2004). 

[11] F.W.J. Olver, Asymptotics and special functions, Academic Press, San Diego 
(1974). 

[12] S. Ponnusamy, M. Vuorinen, Asymptotic expansions and inequalities for hyper- 
geometric functions, Mathematika 44, 278-301 (1997). 

[13] J. Poschel, Quasi-periodic solutions for a nonlinear wave equation, Comment. 
Math. Helvetici 71, 269-296 (1996). 

[14] G.N.Watson, A treatise on the theory of Bessel functions, Reprint of the Sec- 
ond (1944) Edition, Cambridge Mathematical Library, Cambridge Univ. Press 
(1995). 



44 



